High-Precision Relative Locations of Two Microearthquake Clusters in Southeastern Sicily, Italy

In November 1999 and in January 2000, two microearthquake swarms occurred in Southeastern Sicily (Italy). They were analytically located in the depth range 17-25 km, some kilometers northward from the buried front of a regional foredeep, below the active thrust zone of the Sicily mountain chain. Their hypocentral distribution showed two distinct clusters, and comparison of the waveforms revealed clearly that the two swarms formed two distinct families of multiplet events. This led us: i) to carry out a precise relocation relative to two chosen master events of the families, and ii) to better define the geometrical structure of the two clusters. The cross-spectral method was applied to obtain precise readings of the wave onsets. SH-wave onsets were used instead of P-waves, as they showed clearer onsets and a good signal-to-noise ratio. Residuals of the relative locations showed small values, no more than several meters on average. The vertical extent of the two relocated clusters was 500 m and 250 m, respectively, while the horizontal extent was 250 m. Hypocenters of the first cluster clearly delineate a NNW trending plane with almost vertical dip, matching one nodal plane of the focal mechanism obtained as a composite solution of all events of the cluster. Given the considerable gap angles, because of unfavorable network geometry with respect to the events, the stability of our results was tested carrying out a Montecarlo experiment. Varying the onset times randomly in the range of 5 ms, a dispersion of the locations less than 10 m in longitude, and less than 50 m both in latitude and depth was found. Similar results were obtained when comparing relocations carried out with different master events. Thus, the overall geometrical characteristics of the clusters were not affected seriously by random errors. Considering the geo-structural framework of the region, together with the location and time evolution of the two clusters, fluids of plutonic origin are suggested as the trigger mechanism

Seismicity, seismotectonics and crustal velocity structure of the Messina Strait (Italy)

The Messina Strait is the most important structural element interrupting the southernmost part of the Alpine-Apenninic orogenic belt, known as the Calabro-Peloritan Arc. It is being a narrow fan-shaped basin linking the Ionian Sea to the Tyrrhenian Sea. This region is affected by considerable seismic activity which mirrors the geodynamic processes due to the convergence between the African and the Eurasian plates. In the last four centuries, a significant number of disastrous earthquakes originated along the Arc. Among these, the most noteworthy event occurred on December 28, 1908 (known as the Reggio Calabria-Messina earthquake), in the Messina Strait area and caused a large tsunami and more than 100,000 casualties. In this research we focus on the relationships between the general tectonic setting, which characterize the Messina Strait and adjacent areas, seismicity patterns and the crustal structure. We analyzed a data set consisting of more than 300 events occurring in the years from 1999 to 2007, having a magnitude range from 1.0 to 3.8. This data set was exploited in a local earthquake tomography, by carrying out a simultaneous inversion of both the three-dimensional velocity structure and the distribution of seismic foci. We applied the “tomoADD” algorithm, which uses a combination of absolute and differential arrival times and a concept of self-adapting grid geometry, accounting for ray density encountered across the volume. With this method the accuracy of event locations is improved and velocity structure near the source region is resolved in more detail than standard tomography. Fault plane solutions were obtained for the major and best-recorded earthquakes. The obtained velocity images highlight vertical and lateral heterogeneities that can be associated with structural features striking from NNE-SSW to NE-SW. These results are consistent with important tectonic elements visible at the surface and the pattern delineated by earthquake locations and focal mechanisms

Estimation of an optimum velocity model in the Calabro-Peloritan mountains – Assessment of the variance of model parameters and variability of earthquake locations

Accurate earthquake locations are of primary importance when studying the seismicity of a given area, they allow important inferences on the ongoing seismo-tectonics. Both, for standard, as well as for earthquake relative location techniques, the velocity parameters are kept fixed to a-priori values, that are assumed to be correct, and the observed traveltime residuals are minimised by adjusting the hypocentral parameters. However, the use of an unsuitable velocity model, can introduce systematic errors in the hypocentre location. Precise hypocentre locations and error estimate, therefore, require the simultaneous solution of both velocity and hypocentral parameters. We perform a simultaneous inversion of both the velocity structure and the hypocentre location in NE-Sicily and SW-Calabria (Italy). Since the density of the network is not sufficient for the identification of the 3D structure with a resolution of interest here, we restrict ourselves to a 1D inversion using the well-known code VELEST. A main goal of the paper is the analysis of the stability of the inverted model parameters. For this purpose we carry out a series of tests concerning the initial guesses of the velocity structure and locations used in the inversion. We further assess the uncertainties which originate from the finiteness of the available datasets carrying out resampling experiments. From these tests we conclude that the data catalogue is sufficient to constrain the inversion. We note that the uncertainties of the inverted velocities increases with depth. On the other hand the inverted velocity structure depends decisively on the initial guess as they tend to maintain the overall shape of the starting model. In order to derive an improved starting model we derive a guess for the probable depth of the MOHO. For this purpose we exploit considerations of the depth distribution of earthquake foci and of the shear strength of rock depending on its rheological behaviour at depth. In a second step we derived a smooth starting model and repeated the inversion. Strong discontinuities tend to attract hypocentre locations which may introduce biases to the earthquake location. Using the smooth starting model we obtained again a rather smooth model as final solution which gave the best travel-time residuals among all models discussed in this paper. This poses severe questions as to the significance of velocity discontinuities inferred from rather vague a-priori information. Besides this, the use of those smooth models widely avoids the problems of hypocentre locations being affected by sudden velocity jumps, an effect which can be extremely disturbing in relative location procedures. The differences of the velocity structure obtained with different starting models is larger than those encountered during the bootstrap test. This underscores the importance of the choice of the initial guess. Fortunately the effects of the uncertainties discussed here on the final locations turned out as limited, i. e., less than 1 km for the horizontal coordinates and less than 2 km for the depth

Shear wave splitting changes associated with the 2001 volcanic eruption on Mt. Etna

The time delays and polarizations of shear wave splitting above small earthquakes show variations before the 2001 July 17–August 9 2001 flank eruption on Mt Etna, Sicily. Normalized time delays, measured by singular value decomposition, show a systematic increase starting several days before the onset of the eruption. On several occasions before the eruption, the polarization directions of the shear waves at Station MNT, closest to the eruption, show 90◦- flips where the faster and slower split shear waves exchange polarizations. The last 90◦-flip being 5 days before the onset of the eruption. The time delays also exhibit a sudden decrease shortly before the start of the eruption suggesting the possible occurrence of a ‘relaxation’ phenomena, due to crack coalescence. This behaviour has many similarities to that observed before a number of earthquakes elsewhere

Shear wave splitting time variation by stress-induced magma uprising at Mount Etna volcano

Shear wave splitting exhibits clear time variations before the July 17th – August 9th, 2001 flanK eruption at Mount Etna. The normalized time delays, Tn, detected through an orthogonal transformation of singular value decomposition, exhibit a clearincrease starting 20 days before the occurrence of the eruption (July 17th); the qS1 polarization direction, obtained using a 3D covariance matrix decomposition, shows a 90°-flip several times during the analyzed period: the last flip 5 days before the occurrence of the eruption. Both splitting parameters also exhibit a relaxation phase shortly before the starting of the eruption. Our observations seem in agreement with Anisotropic Poro Elasticity (APE) modelling, suggesting a tool for the temporal monitoring of the build up of the stress leading to the occurrence of the 2001 eruption at Mt. Etna

PRIN Project 2010-11 “Active and recent geodynamics of Calabrian Arc and accretionary complex in the Ionian Sea”: new constraints from geological, geodetic and seismological data

This contribution illustrates the preliminary results of our Research Unit in the PRIN Project 2010-11, which focuses on active and recent geodynamics of Calabrian Arc. The integration of the new geological, geodetic and seismological data supports the inferred recent plate boundary reorganization in the central-southern Mediterranean, where the regional GNNS velocity fields point to a deceleration or cessation of Calabrian Arc migration, and to extension along the axis of the Calabrian Arc, accommodated by normal faulting (e.g. Capo Vaticano and Messina Straits (Aloisi et al., 2012; Pepe et al., 2014; Spampinato et al., 2014). The study of the lateral borders of the Arc revealed that oblique strike-slip displacement has occurred during its southeastwards migration. Active dextral transtension is occurring along the NNW-striking Aeolian-Tindari Letojanni fault system, forming the southern boundary of the Arc. It joins to the north other two boundaries characterized by different tectonic regimes, a contractional belt in the southern Tyrrhenian sea, where a tectonic inversion has occurred since the middle Pleistocene, and the extensional one in northeastern Sicily and western Calabria (Palano et al., 2012; Barreca et al., 2014a). Along the northern boundary of the Arc, the so-called Pollino line (onshore) and Sibari Line (offshore), active deformation has been documented on folds growing above blind oblique thrust ramps extending offshore, controlling the present morphobathymetric pattern (Santoro et al., 2013). Although external to the Calabrian Arc, we also devoted attention to the front of the Maghrebian thrust belt in western Sicily where we presented the first evidence of historical co-seismic deformation on a thrust array running from the Belice area to the Sicily Channel (Barreca et al., 2014b). Morphotectonic analysis and fault numeric modeling of uplifted Pleistocene marine terraces and Holocene paleo-shorelines has documented that most of the uplift along the Calabrian Arc is related to regional processes and the residual to coseismic displacement on major faults, both transpressional and transtensional, at the borders, and extensional along the chain axis

Accurate hypocentre locations in the Middle-Durance Fault Zone, South-Eastern France

A one-dimensional velocity model and station corrections for the Middle-Durance fault zone (south-eastern France) was computed by inverting P-wave arrival times recorded on a local seismic network of 8 stations. A total of 93 local events with a minimum number of 6 P-phases, RMS<0.4 s and a maximum gap of 220° were selected. Comparison with previous earthquake locations shows an improvement for the relocated earthquakes. Tests were carried out to verify the robustness of inversion results in order to corroborate the conclusions drawn from our findings. The obtained minimum 1-D velocity model can be used to improve routine earthquake locations and represents a further step toward more detailed seismotectonic studies in this area of south-eastern France

Multidisciplinary study of the Tindari Fault (Sicily, Italy) separating ongoing contractional and extensional compartments along the active Africa–Eurasia convergent boundary

The Africa–Eurasia convergence in Sicily and southern Calabria is currently expressed by two different tectonic and geodynamic domains: thewestern region, governed by a roughlyN–S compression generated by a continental collision; the eastern one, controlled by a NW–SE extension related to the south-east-directed expansion of the Calabro–Peloritan Arc. The different deformation pattern of these two domains is accommodated by a right-lateral shear zone (Aeolian–Tindari–Letojanni fault system) which, from the Ionian Sea, north of Mt. Etna, extends across the Peloritani chain to the Aeolian Islands. In this work, we study the evidence of active tectonics characterizing this shear zone, through the analysis of seismic and geodetic data acquired by the INGV networks in the last 15 years. The study is completed by structural and morphological surveys carried out between Capo Tindari and the watershed of the chain. The results allowed defining a clear structural picture depicting the tectonic interferences between the two different geodynamic domains. The results indicate that, besides the regional ~N130°E horizontal extensional stress field, another one, NE–SW-oriented, is active in the investigated area. Both tension axes are mutually independent and have been active up to the present at different times. The coexistence of these different active horizontal extensions is the result of complex interactions between several induced stresses: 1) the regional extension (NW–SE) related to the slab rollback and back-arc extension; 2) the strong uplift of the chain; 3) the accommodation between compressional and extensional tectonic regimes along the Aeolian– Tindari–Letojanni faults, through a SSE–NNW right-lateral transtensional displacement. In these conditions, the greater and recurring uplift activity is not able to induce a radial extensional dynamics, but, under the “directing” action of the shear system, it can only act on the regional extension (NW–SE) and produce the second system of extension (NE–SW)

Assessing seismic efficiency from scalar Moment-rates: an application to Mt. Etna volcano (Italy)

Here we propose an improved estimation of the scalar seismic (from instrumental and historical catalogues), geodetic and geologic moment-rates for the eastern flank of Mt. Etna. The estimated moment-rates have been compared in terms of seismic efficiency. Results show that all the calculated efficiency values are lower than 40%, i.e., the geodetic moment-rate estimations are generally larger than the seismic and the geologic ones. Although a number of reasons may account for the observed discrepancy, we are confident that a large amount of the deformation affecting the eastern flank occurs aseismically

Terremoti con effetti macrosismici in Sicilia orientale nel periodo Gennaio 2009 - Dicembre 2013

Sono analizzati ed elaborati i dati dei rilievi macrosismici relativi ai terremoti verificatisi in Sicilia nel periodo 2009-2013. Si tratta di eventi che hanno provocato effetti macrosismici di rilievo e/o danneggiamento per la maggior parte localizzati nellâ area etnea, cui si aggiungono alcuni terremoti di magnitudo moderata, legati a sequenze sismiche significative verificatisi nei settori ibleo, peloritano e eoliano. I dati sono stati raccolti e elaborati secondo le procedure operative di prassi adottate in questi casi dal gruppo QUEST (ex-TTC 1.11 â Osservazioni e monitoraggio macrosismico del territorio nazionaleâ ) dellâ INGV, e successivamente parametrizzati secondo gli standard adottati per la compilazione del catalogo CPTI e banca dati macrosismica DBMI (ex-TTC 5.1 â Banche dati e metodi macrosismiciâ )